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Renormalization factors are most easily extracted by going to the massless limit of the quantum 
field theory and retaining only a single momentum scale. We derive the factors and renormalized 
Green functions to all orders in perturbation theory for rainbow graphs and vertex (or scattering) 
diagrams at zero momentum transfer, in the context of dimensional regularization, and we prove 
that the correct anomalous dimensions for those processes emerge in the limit D —* 4. 
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. I. INTRODUCTION 

t> , 

The connection between knot theory and renormaliza- 
tion theory Q is one of the more exciting developments 
of field theory in recent years because it relates appar- 
ently different Feynman diagrams through the common 
topology of the associated knots. Thus it serves to ex- 
plain why transcendental numbers for the renormaliza- 
tion constants Z occur in some diagrams Q and not 
in others, thereby allowing the Feynman graphs to be 
grouped into equivalence classes. Sometimes, in gauge 
theories the Z factors within a particular class may can- 
cel because of the existence of Ward identities, leaving a 
non-transcendental result for Z; this happens in electro- 
dynamics of scalar and spinor particles to fourth order 
in the quenched limit and in chromodynamics to third 
order ||. 

The class of graphs which correspond to ladders and 
rainbows are especially simple in this connection because 
they possess trivial knot topologies. Thus one may an- 
ticipate that Z-factors for them are particularly easy to 
evaluate. Kreimer Q| has provided rules for extract- 
ing them within the framework of dimensional regular- 
ization, through the standard expedient of finding the 
simple 1/e = 2/(4 — D) pole term arising in products 
of functions, after removing lower-order pole terms con- 
nected with subdivergences. Thus vertex diagrams bring 
in function factors of the type 

, A(e) = (p 2 )^ +1 > 



while rainbow graphs lead to products of 
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It has to be said that, although the procedure is straight- 
forward, extracting the 1/e term in n'th order requires 
considerable graft. Kreimer has proven that the simple 
pole in e is free of Riemann zeta-functions. 

In this paper we shall show that the problem can be 
solved to all orders in perturbation theory for ladders 
and rainbows Q, in the context of renormalization in 
dimensional regularization because of two fortuitous cir- 
cumstances: (i) the Green function satisfies a differen- 
tial equation and (ii) this equation is actually soluble 
in terms of Bessel functions. The limit as D — » 4 may 
then be taken at the end and, as a useful check, the 
anomalous dimension properly emerges. (It is a rather 
delicate limit, requiring a saddle-point analysis of the 
integral representation of the Bessel function, since it 
looks quite singular.) We have successfully carried out 
this program for meson- fermion theories, both for vertex 
functions and rainbow diagrams; however we have not 
succeeded in solving the problem near D — 6 for (f> 3 the- 
ory because the differential equation is of fourth order 
and cannot be expressed in terms of standard functions; 
nevertheless we can obtain the answer in the limit x — > 
or p — > oo for D — 6. 

In the next section we treat the vertex diagrams for 
scalar mesons, while the following section contains the 
analysis of the rainbow diagrams. The appendix contains 
details of the vector meson case, which are rather more 
complicated. 



II. VERTEX DIAGRAMS 



We shall consider a theory of massless fermions and 
mesons <p in D-dimensions since the purpose of our work 
is to investigate the behaviour of the Green function 
as D tends to 4. Let 7[ r ] signify the product of r 7- 
matrices (of size 2 D / 2 x 2 D/>2 ) normalized to unity, namely 
7[^i^ 2 - /Jr] S0 that we can write the meson- fermion inter- 
action in the form 

Ant = gtpj[ r ]1p(j) [r] , 

where is the corresponding tensor meson field. The 
equation for the renormalized tensor vertex function Pui 
at zero meson momentum, taking out the factor <?, is 

r W (p) = Zj [r] S r s - ig 2 J d D q^ [r] -^-T [s] (q)-^--f [r , ] A rr (p-q). 

(1) 

We shall assume that the massless meson propagator 
above, A rr (p — q), can be chosen in a Fermi- Feynman 
gauge so 

A"-' = (-l)V r 7(p-9) 2 , 
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where r/ stands for the diagonal Minkowskian met- 
ric pertaining to the tensor structure, specifically 

To make further progress we utilize the non-amputated 
Green function, 

1 1 

G[ r ](p) = — r M 0) — , 

J.p "f.p 
to remain with the 'simpler' linear integral equation, 

1 .pG [s] {p) 1 .p = Z5 r sl[r] + i(-Yc s r g 2 J d D g (2) 

where dPq = d D q/ (2ir) D . The nature of the couplings in 
massless theories means that the Green function always 
stays proportional to 7[ s ] and can be decomposed into 
just two pieces H, 

G[ s ](p) = i[ s ]A{p 2 ) + -f.pj ls] j.pB(p 2 ). (3) 

On the right-hand side of (2), c*7f s ] = 7[ r ]7[ s ]7^, is es- 
sentially an element of the Fierz transformation matrix 
for any D, given by || 

c; = (-ir X>i)«(^ 9 r )Q. 

We shall convert the integral equation (2) into a differ- 
ential equation by taking the Fourier transform. (In fact 
we could almost have done this from the word go by writ- 
ing the equation for the full Green function in coordinate 
space.) This manoeuvre produces 

1 .dG [s] {x) 1 ^d = ZS r sl[r] 5 D (x)+ic s r (-) r g 2 A c (x)G [s] (x), 

(4) 

where the massless meson propagator is iA c (x) = 
T{D/2~l)(-x 2 + ie) 1 - D / 2 /4TT D / 2 . Because the coupling 
constant is a dimensionful quantity, we can define a di- 
mensionless strength 'a = fine structure constant /47r', 
via 

(-) r c s r g 2 T{D/2 - 1) = 16Tr D / 2 ^- D a 

upon introducing a mass scale fi. This simplifies the re- 
sulting expressions, as we can see in the purely scalar 
case, where there is but a single term and equation: 

G{x) = -Z5 D {x). (5) 

The equation is readily solved in dimension D = 4 yield- 
ing G CX (-X 2 +ig)(-l-Vl+4a)/2 _ For D 4 we can make 

progress by passing to a Euclidean metric (r 2 = —x 2 ), 

G{r) = ZS D {r). 



Q2 _^_ { _ V 2 X 1 ) 2-DI2 



dr 2 



D-l d 
r dr 



4a 4 _ D 
^2 ^ r > 
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The whole point of the manipulation is that one is for- 
tunately able to solve this equation (for r ^ at first) in 
terms of known functions, namely Bessel functions. The 
correct choice of solution, up to an overall factor, is 



G(r) oc r-V^^V^tHVe); D = 4-2e. 

because in the limit as a — > we recover the free-field 
solution r 2 ~ D . For dimensional reasons, let us carry out 
our renormalization so that G(l//i) = /i 2 ~ 2£ . With that 
convention, the vertex function reduces to 

G(r) = (r/^Y~ 1 J 1 _ 1/e (V^(firY/e)/J 1 _ 1/e (V^/e). 

(6) 

Furthermore the correct singularity for the time-ordered 
function S(x) emerges if we reinterpret r 2 = —x 2 + ie 
above. We shall not worry at this stage whether a is 
positive or negative — the sign can vary with the model 
anyway — since we can easily continue the function from 
J to I as needed. 

The problem presents itself: how does the four- 
dimensional result, with its anomalous scale 7 = 
VI + 4a — 1, emerge from (6) as e — > 0~, say? This 
is clearly a delicate limit because both the index and the 
argument of the Bessel function become infinitely large. 
Before answering this question, let us note that in a per- 
turbative expansion of (4), viz. a small argument expan- 
sion of J, 



[1 

G(r) = r 2e - 2 



q(^) 2g , a 2 (^) 4 ' 

c(2e-l) 2e 2 (2e-l)(3e-l) 



[1 + e(2 e a -l) + 2£ a (2e-l)(3e-l) + " "] 

the poles in e cancel out to any particular order in a. For 
instance up to order a 2 we obtain as e — > 0, 

G(r) = r- 2 [l + 2(-a + a 2 ) In fir + 2a 2 (ln^r) 2 + • • •] 

which agrees precisely with the expansion of the anoma- 
lous dimension in the logarithmic terms. Returning to 
the limit of small e, we will make use of the saddle point 
method of obtaining asymptotic expansions of integrals. 
Suppose that f(t) has a minimum at t = r in the integral 
representation, 



f = it f ex P 



Then the saddle point method gives 



F = cxp/(r) ' 



a/27t/"(t) 



1+ f""{r) 5(/"'(r)) 2 , 35(/""(r)) 2 



8(/"(r)) 2 24(/"(r))3 384(/"(r)) 4 



As confirmation of the correctness of the terms above 
we can verify that the Debye expansion |)| of the Bessel 
function is properly reproduced, 



J„(tVcoshr) = ^ J 



27T + 00 

v(sinhi/ coshr — t) 

coth r ( 5 , , n 
1 + 1 - -(cothr) 2 



V27r^tanhr 



8^ V 3 
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because the integrand minimum occurs at t = r. Our 
case (6) is a variant of this. Working only to first order in 
e, and taking a negative initially, the integrand exponent 
is 

(V— 4a sinhi + t)/e + \J — 4a In fir sinhi — < 
and is stationary at the complex value t = r, where 

cosh r 



V-4a(l + eln/xr) 

Following through the mathematical steps, and omitting 
straightforward details, to order e we end up with, 

2ae , , 

1 -TT^ ln(H + '-- 

(7) 



G(r) ^r^ifxr) 1 



It is satisfying that this produces the all-orders (in cou- 
pling, a) result at 4-dimensions when e — > 0, with the 
correct anomalous scale. 

The problem can be treated in much the same way 
for a pseudoscalar meson field. The only possible dif- 
ference is a change in sign of a, because of '75' matrix 
anticommutation. As for the vector case (r = s = 1), 
the procedure produces a pair of coupled equations for 
the two scalar components A and B of the Green func- 
tion, G^{p) — j^A(p 2 ) + j.pj fl j.pB(p 2 ). A discussion of 
this case is given in the appendix, where it is shown that 
the only easy limit is D = 4; one finds after Euclidean 
rotation that 

A = ar 13 , B = br< 3+2 ; a/b = c(f3 + 2)/{p - 2), 

where 



= -1 + V 5 + V16 + 4c + c 2 and c = g 2 /2ir 2 , 

in four dimensions; we have chosen the root which re- 
duces to the free field solution when g = although one 
can contemplate strictly non-perturbative solutions [fl0|| . 

The difficulty is symptomatic of what happens in <fi 3 
theory near six dimensions; in that case the Green func- 
tion, G(p) — T(p)/p i obeys the Fourier transformed 
equation |j, 



<9 4 -4a 



(pr) 



6-D 



G(x) = ZS D (r). (8) 



This is a differential equation of fourth order in r and its 
solution cannot readily be expressed in terms of familiar 
transcendental functions. However in the limit as D — > 6, 
it is quite simple to find the (power law) solution: 



G(r) oc r''; (3 = -1 - y5 + 2V4+~a, 

which correctly reduces to the free-field solution (3 = — 4 
when the coupling vanishes. 
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III. RAINBOW DIAGRAMS 



We wish to treat the corrections to the fermion prop- 
agator in a similar manner, by considering the rain- 
bow corrections. In such an approximation the rainbow 
graphs give rise to a self-energy, which is self-consistently 
determined according to 



£r(p) = -ig 2 J '■ 



d D q 



1 1 v ( n \ 1 
£fl(£?) — 

7-9 7-9 7-9. 



(p - q) 2 

with the unrenormalized propagator determined by 

Sr(p) = — — £r(p) — 

at this level [B. This leads to the renormalized rainbow 
corrected propagator equation, 

l-pS R (p)j-P= Z^-f.p + ig 2 J j—^yS(q). (9) 

The nature of the massless problem is that one can always 
write Sr(p) — j.pa(jp), and by Fourier transformation, 
convert (9) from an integral equation to a differential 
equation, 

-ij.dd 2 a{x) = i-i.dZ^5 D {x) + ig 2 A c {x) ij.da(x) 

or 

[d 2 + ig 2 A c (x)] d^a(x) = -Z^S D (x) (10) 
Now for any function /(Vx 2 ), using the two lemmas, 

d„f = xJ'/Vx*, (11a) 

W = U v - ^) -£= + ^/", (lib) 



we can carry out an Euclidean rotation in order to arrive 
at the differential equation for the scalar function S = 
da/dr: 

d{D-l,d\ , 4a(H— s{r)= _ z ^ 

(12) 



dr V r dr 



In 4-D this has the simple solution S(r) oc r 1 2v/1+a , 
in turn implying S(p) oc 7.pp~ 4+2%/1+a . However it is in 
fact possible to solve (12) for any D. The proper solution, 
normalized to 5(l//i) = /i 3_2e is 



_ 2 J 1 _2/ e (V = S(H7e) 



5(r) - r(r/^- 2 ^/ ' '- . (13) 

Ji-2/e(v-4a/e) 

To get the rainbow propagator, we must first integrate, 
= fS{r)dr oc ELo^-^+smlV^H'/e). 
and then Fourier transform to obtain Snip) = ^.pa(p). 

If the mesons are neither scalar nor pseudoscalar, but 
tensor, the coupling constant is multiplied by the factor 
c\; that is all. 
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IV. CONCLUSIONS 



We have demonstrated that it is possible to work out 
the all-orders solution of Green functions for ladder and 
rainbow diagrams for any dimension D and that, in the 
limit as D approaches the physical dimension, the correct 
scaling dimension is obtained. We have exhibited fully 
how this happens for scalar theories, but have succeeded 
only to a limited extent in vector theories, because the 
equations are coupled and end up as fourth order ones, 
with no transparent expression in terms of standard func- 
tions of mathematical physics. In any event, it is clear 
from the form of the Green function that there are no 
transcendental constants in sight, even when we expand 
the answers perturbatively in terms of ln(/ir), so that the 
renormalization constants are free of them. This confirms 
the finding of Kreimer for arbitrary ladder/rainbow order 
H and does not come as a surprise. 

One can extend the ideas here to scattering processes 
which contain a single momentum scale, such as fermion- 
fermion scattering (again ladder graphs) for any D. It is 
a simple matter of taking the Fourier transform in par- 
ticular channels and converting the momentum integral 
equations to differential ones in coordinate space. We 
shall not labour the issue in this paper since the steps are 
fairly obvious and can easily be filled in by the reader. 
What we have not solved for any D is the case of crossed 
ladders, when the kernel will presumably lead to tran- 
scendental Z constants; that is a task for the future. 
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APPENDIX: THE VECTOR CASE 



The vector vertex function 

upon Fourier transformation and tracing with 7„, pro- 
duces the coordinate space equation, 

(d^^d^A+d^B = Zr 1 ^5{x)+ig 2 {D-2)A c {x)[r 1 ^A+{d\ liv -2d ll d v )B}. 

Using lemmas (11), and identifying the terms multiplying 
rj^ v and x^x v , we arrive at the pair of coupled equations, 



o a A 

0[-A + OB}+^ = zg 2 (D-2)A c 
2 ar 



' 2dB 
A+- — 

r ar 



OB 



QA = ig 2 (D~2)A c QB, 



where O = ( 



p-i d 

dr ' 



and Q = ( 



1 d_ 

r dr I 



We 



have not suceeded in solving these equations in terms of 
familiar functions for D ^ 4. However in 4-D, one can 
make considerable progress by looking for a power-law 
solution of the type, A{r) = ar 13 , B(r) = br@ +2 . Simple 
calculation reveals that a solution exists provided that 

a/b = c{0+2)/{0-2), c 2 +4c-(3(f3-2){f3+2)(f3+4) = 0, 

The quartic in the power exponent (3 is fortunately simple 
to solve in terms of the coupling (or c), the answer being 



(3 = -l+ v 5 + V16 + 4c + c 2 = 2 + c/3 - c 2 /54 • ■ 
so that a/b~ 12 + 5c/3 H . 



g 2 /2n 2 
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